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Optimal conversion of Bose-Einstein condensed atoms into molecules via a Feshbach 
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In many experiments involving conversion of quantum degenerate atomic gases into molecular 
dimers via a Feshbach resonance, an external magnetic field is linearly swept from above the res- 
onance to below resonance. In the adiabatic limit, the fraction of atoms converted into molecules 
is independent of the functional form of the sweep and is predicted to be 100%. However, for non- 
adiabatic sweeps through resonance, Landau-Zener theory predicts that a linear sweep will result 
in a negligible production of molecules. Here we employ a genetic algorithm to determine the func- 
tional time dependence of the magnetic field that produces the maximum number of molecules for 
sweep times that are comparable to the period of resonant atom-molecule oscillations, 27ffi^ aM . The 
optimal sweep through resonance indicates that more than 95% of the atoms can be converted into 
molecules for sweep times as short as 27rn^ bi while the linear sweep results in a conversion of only 
a few percent. We also find that the qualitative form of the optimal sweep is independent of the 
strength of the two-body interactions between atoms and molecules and the width of the resonance. 

PACS numbers: 03.75.-b,03.75.Pp 
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I. INTRODUCTION 

Unlike alkali atoms, molecules can not be directly 
cooled using the laser cooling techniques that led to Bose- 
Einstein condensation because of the complex rotational- 
vibrational spectrum of the molecules. As a result, 
two-photon Raman photoassociation and Feshbach res- 
onances have become the standard tools to create trans- 
lationally cold molecules starting from ultra-cold atomic 
gases. The conversion of a macroscopic number of quan- 
tum degenerate atoms into molecular dimers starting 
from either a Bose-Einstein condensate |U- LJ Lj- Q or a 
Fermi gas [a S 13 nas been observed by several exper- 
imental groups using a Feshbach resonance. This work 
culminated in the formation of a molecular Bose-Einstein 
condensate (MBEC) HG3- 

Until recently the field of molecular optics was in the 
same state as atom optics before BEC. Diffraction and 
interferometry had been demonstrated lllj but there was 
no source of high density phase coherent monoenergetic 
molecules analogous to a laser that could be used to ob- 
serve nonlinear and quantum optical effects. Indeed, the 
ability to now coherently produce molecules with a high 
phase space density opens up new avenues of research in 
the area of matter- wave optics such as lasing and matter 
wave amplifications with molecular fields, nonlinear mix- 
ing of atomic and molecular matter waves, and the gen- 
eration of entangled atoms by 'down conversion' of the 
molecules. Recent experiment have shown the phase co- 
herent and momentum conserving nature of matter wave 
second harmonic generation 12]. 

From the perspective of atom optics, the conversion 
of atoms into molecules via a Feshbach resonance or 
photo-association is the matter-wave analog of second 
harmonic generation of photons in a nonlinear crystal 
with a x^ susceptibility, which has been used to cre- 



ate entangled photon states. From another perspective, 
these methods can be viewed as the first step towards 
'quantum super-chemistry' where chemical reactions oc- 
cur in a controllable phase coherent manner and exhibit 
novel quantum features such as interference and bosonic 
stimulation|l3|. Also, the ability to create molecules with 
permanent dipole moments using a heteronuclear Fesh- 
bach resonance [14j opens up the possibility of studying 
bosonic systems with anisotropic interactions 15]. Dipo- 
lar molecules can be stored in electrostatic storage rings 
[16j that have an area that is twenty five times larger 
than the largest neutral atom magnetic storage ring [ljj • 
Storage rings for dipolar molecules therefore have the 
potential to be used for high precision rotation sensors 
based on the Sagnac effect, which is proportional to area 
enclosed by the ring. 

Whatever the reason for creating ultracold molecules, 
it is highly desirable to have an efficient means of pro- 
duction. While the molecules created via a Feshbach 
resonance are translationally very cold, they are vibra- 
tionally very hot and can decay to lower lying vibrational 
states via exoergic inelastic collisions with atoms or other 
molecules. For an atomic BEC, the molecular two-body 
decay rates are of the order 1CP 11 — 10~ 10 cm 3 /s, which 
gives a lifetime of 100/is for a typical atomic density 
0, EH Il3 • It is therefore important to be able to create 
the molecules as quickly as possible so that experiments 
can be performed on them before they are lost from the 
system. (In this paper, we will not consider molecules 
created from an atomic Fermi gas, which can have life- 
times on the order of a second due to Pauli blocking of 
collisions |20j.1 

In experiments using Feshbach resonances in atomic 
BECs, adiabatic rapid passage is used to convert atoms 
into molecules by 'slowly' varying an external magnetic 
field from above resonance, where the energy of the 



molecular state lies above that of the atoms to below 
resonance, where molecules are the stable ground state. 
When the magnetic field strength is swept linearly across 
the Feshbach resonance, the atomic and molecular pop- 
ulations can be predicted by the Landau-Zener (LZ) for- 
mulae for a two-level quantum system 21] . If there 
are initially no molecules at t — » — oo, the probability 
at t — > +00 of a pair of atoms remaining unconverted, 
Patom, and of being converted into molecules, P mo iecule, 
is given by |21j. 



1 atoms 6 

= 1 



-2k&, 



molecules 



'LZ 



-2tt8 l 



«L«/4|A| 



(1) 

(2) 
(3) 



where flR a bi is the Rabi frequency coupling the two states 
and A(i) = |A|i is the energy difference (detuning) be- 
tween the two states. The Landau-Zener (LZ) factor Slz 
is the main parameter that describes the atom-molecule 
conversion in theoretical models 22] . 

In the absence of decay of the atoms and molecules, 
100% molecular conversion is expected from the LZ the- 
ory for sufficiently small | A| . If the detuning varies slowly 
enough to satisfy the adiabatic condition, 5lz > 1, the 
initial atomic BEC will adiabatically evolve into a molec- 
ular BEC. In reality, the conversation rates are always 
much smaller than that predicted by LZ theory due to 
molecular vibrational decay \loL [lj| • Experiments with 
atomic BECs have yielded molecular conversion efficien- 
cies of only ~ 5 — 10% with a large fraction of the atoms 
being lost during the sweep [3, 3 • The missing fraction 
of atoms in the experiments are attributed to vibrational 
decay of the molecules, although spontaneous dissocia- 
tion of molecules caused by inelastic spin flips has been 
shown to be important in 85 Rb 24]. This was recently 
confirmed by experiments in an optical lattice where 87 i?6 
dimers were created with 95% efficiency on lattice sites 
initially containing only two atoms so that vibrational 
decay of the resulting molecule was impossible |25j . 

Motivated by the need to create molecules faster than 
atom and molecule loss rates, we have explored atom- 
molecule conversion via a Feshbach resonance in the non- 
adiabatic regime. The sweep times that we are primarily 
interested in are shorter than the molecular lifetimes and 
therefore we ignore decay processes. To this end we have 
employed a genetic algorithm to determine the magnetic 



field sweep that maximizes the fraction of atoms con- 
verted into molecules for sweeps as short as 27rf2^ bi . 
We find that the optimal sweep shows a conversion ef- 
ficiency in excess of 95% when the linear sweep results in 
a conversion of only a few percent. Moreover, the general 
form of the optimal sweep is independent of the strength 
of the resonance and the strength of the two-body in- 
teractions between atoms and molecules. These results 
indicate that higher conversion efficiencies could be ob- 
tained in experiments with the added bonus of longer 
times to perform experiments on the molecules. A recent 
experiment on a Cs atomic BEC showed that up to 30% 
of the atoms could be converted into molecules using a 
sweep through resonance that is similar to the one we 
propose here [26| . However, the duration of their sweep 
is much longer than what we consider here and can be 
considered adiabatic. 

This paper is organized as follows. In section II we 
explain our physical model and in section III we discuss 
the genetic algorithm used. In section IV we compare 
the results obtained from the genetic algorithm with a 
linear LZ sweep through resonance. Section V presents 
some conclusions and future directions for this work. 



II. PHYSICAL MODEL 

A Feshbach resonance is a collisional resonance be- 
tween a pair of free colliding atoms and a molecular 
bound state of those atoms. The coupling between the 
atomic and molecular states is due to the hyperfine in- 
teraction, which can result in spin flip of the electrons 
|27|. The colliding atom state is usually referred to as 
the open channel while the molecular state is known as 
the closed channel since it can only be accessed by a spin 
flip of one of the atoms. Because the atomic and molec- 
ular states have different magnetic moments, the energy 
difference of the atoms and moleculs can be tuned using 
an external magnetic field. 

At zero temperature, a weakly interacting BEC of 
atoms can be described by the semi-classical Gross- 
Pitaevskii equation. These equations can readily be ex- 
tended to include a molecular BEC produced by coherent 
interconversion of atoms into molecules via a Feshbach 
resonance El El HI, 
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where (f>i is the atomic BEC wave function and </> 2 that of the molecules. Here, to* are the masses and Vi(x) 



the external trapping potential of the atoms (i = 1) and 
molecules (i = 2). The two-body interaction is given by 
a contact potential, V(r — r') = Uij8(y — r'), with the 
coupling constants Uij — 2-kTi dij/niij where an is the 
atom-atom s-wave scattering length, 0,22 the molecule- 
molecule scattering length, a\ 2 the atom-molecule scat- 
tering length, and m.^ = rriimj j '(rrii + rrij)- X IS the 
coupling between the open channel and the closed chan- 
nel defined by x = [27r?i 2 aiiA/iAS/m] 1 / 2 , where A/i is 
difference between magnetic moments of an atomic pair 
and a molecule and AB is the width of resonance 27] . 

HA(t) = Afi(B(t) — Bo) is the Zeeman energy dif- 
ference between the molecules and atoms such that Bq 
is the magnetic field strength at which the energy of 
the molecule equals the open channel collision thresh- 
old. For A(i) < the atoms are the lowest energy 
state while for A(i) > 0, the molecular state is the 
lowest energy state of the system. Although we ignore 
vibrational decay of the molecules in our calculations, 
it can be easily included by making the substitution 
A(i) -► A(t) - in a \Mt)\ 2 ~ iKm\Mt)\ 2 E3 

Unless stated otherwise we used the following values 
for the two-body interactions, Uu/ft — 4.96x 10~ 17 m 3 /s, 
U 12 /h = -6.44 x f(r 17 TO 3 /s, U 22 /h = 2.48 x l(T 17 m 3 /s 
corresponding to the scattering lengths of an 87 Rb con- 
densate. For the atom-molecule coupling constant we use 
x/Ti = 1.91 x 10~ 5 m 3 / 2 /s in all of our simulations |22| . 
We also take V(x) = 0, which implies that <f>i and cf>2 
will be spatially homogenous. This is justified even for 
trapped gases in the local density approximation 29] . 

In this paper, we are interested in the functional form 
of A(t) that optimizes the conversion of atoms into 
molecules starting from the initial condition </>i = y/na, 
<f>2 = 0, and A(t) < 0. We also impose the constraint 
that the initial, A(-T/2), and final, A(T/2), values of 
the detuning, and the sweep time, T, are fixed. We em- 
ploy a genetic algorithm to find the optimal sweep for 
range of different T. The conversion efficiency of the 
sweep is given by the probability that a pair of atoms 



becomes a molecule, P, 



molecule 



2\MT/2)\ 2 /n . 



III. GENETIC ALGORITHM 

Genetic algorithms (GAs) have become a widely used 
tool for solving optimization problems that depend on a 
large number of variables [30j. These multidimensional 
optimization techniques proceed by parameterizing an 
objective function in terms of a finite set of coefficients, 
the chromosome. The genetic algorithm operates on a 
set of chromosomes, which represent the members of a 
population. 

The goal is to find the chromosome that is the global 
minimum of the objective function. In order to find the 
optimal chromosome, GAs employ Darwin's principle of 
survival of the fittest. At each generation, the chromo- 
somes in the population compete with each other for 
survival. The members of the population that do the 



best job of minimizing the objective function, survive to 
the next generation while those that do a poor job are 
eliminated. The surviving members are allowed to pro- 
duce offspring for the next generation by some prescribed 
mating rules that mix elements from their chromosomes 
(crossover). Random mutations are also introduced into 
the offspring chromosomes at each generation. This pro- 
cess repeats itself until the change in the optimum so- 
lution between successive generations is less than some 
threshold value. 

We want to optimize the molecular fraction obtained 
by sweeping A(i) from A(-T/2) = -|A | to A(T/2) = 
|Ao|. The detuning is parameterized in terms of either a 
power series, 
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or a Fourier series, 
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The coefficients C n or B n are the elements of a chromo- 
some, where we used N = 160 coefficients for Eq. JSJ and 
N = 320 coefficients for Eq. J7J). The initial population 
consists of N pop = 20 randomly chosen chromosomes. We 
solve the coupled Gross-Pitaevskii equations, Eqs. 01 and 
[5J for each detuning function in the current generation. 
The half of the population that produces the largest frac- 
tion of molecules are kept and allowed to breed to pro- 
duce Np p/ 2 new chromosomes while the other half are 
discarded. 

Crossover is accomplished by choosing pairs of chro- 
mosomes and representing them as binary strings. A 
random point in the strings is chosen and all information 
to the right of that point is swapped. The mutation op- 
erator randomly flips the value of single bits within the 
binary strings of the offspring chromosome. This pro- 
cess of selection, crossover, and mutation continues until 
the difference in optimal conversion efficiency between 50 
successive generations is less than 0.1%. This stopping 
condition is because the GA can become stuck for several 
generations at local minima of the objective function. 



IV. NUMERICAL RESULTS 

The conversion efficiencies of atoms into molecules de- 
termined by the GA are characterized by the initial and 
final values of detuning functions (|Ao|), sweep time 
(T), and the mean field energy shifts for the atoms and 
molecules, which are proportional to Un, U\2, and U 2 2- 
In order to make our results independent of the strength 
of the resonance, all times and frequencies are scaled rel- 
ative to the initial Rabi frequency, Tn a bi/2ir — ^l^ abi = 
h/x^/no- However, to provide a feeling for the relevant 
time scales, we note that T^aU — 16/is for the initial 



atomic density of tiq — 4 x 10 14 cm 3 that we use in our 
simulations. 

We first investigated the molecular conversion effi- 
ciency as a function of magnetic field sweep time with 
|Aq| = 50f2 fiahi . This initial value of the detuning is large 
enough to decouple the atoms and molecules. The cir- 
cles and squares in Fig. 1(a) indicate the conversion effi- 
ciency for a linear sweep and optimal sweep, respectively. 
For long enough sweep times, r = T jT^abi > 100, both 
sweeps show very large conversion efficiencies of around 
95% or more. In this limit, the sweeps are adiabatic and 
the conversion efficiency is therefore independent of the 
shape of the detuning function. 

As t is decreased, the optimal sweep is still able to 
convert almost all atoms into molecules while the con- 
version efficiency of a linear sweep decreases rapidly. For 
t ~ 1, the conversion efficiency of a linear sweep is only 
a few percent. This is easily understood by noting that 
as t decreases the linear sweep is no longer adiabatic. 
For r = 100, 5lz = 1-57 and according to Eq. 0, 



P, 



molecule 



0.99 while for r 



1 the LZ factor de- 



creases to 5lz = 0.015 giving a molecular fraction of 
only Pmoiecuie = 0.089. These values agree qualitatively 
with the numerical solution of the GP equations in Fig. 
1(a), which indicates that the LZ factor is the main pa- 
rameter controlling the atom-molecule conversion. By 
contrast, the optimal sweeps show excellent conversion 
even for very small sweep times. For comparison the 
optimal sweep gives 98.48% and 99.82% conversion for 
r = 1 and 2, while the linear sweep gives 6.74% and 
10.64%, respectively. 

Fig. 2 shows the optimal sweep for r = 2 using both 
Eq. © and J7J). The power series and the Fourier series 
result are almost identical except for the small noise fluc- 
tuations present in the Fourier series. The noise in the 
Fourier series GA varied from run to run but always had 
an average amplitude less than £lR a bi- The noise was the 
result of the high frequency terms in the Fourier series, 
nn/T 3> ^Rabi-, |Ao | , which oscillate much faster than all 
other time scales in the Gross-Pitaevskii equations. We 
therefore do not expect these fluctuations to have any ap- 
preciable effect on the dynamics. These high frequency 
terms needed to be included in A(i) to produce the large 
slope at the beginning and end of the sweep. 

To study the effect of these fluctuations, we superim- 
posed randomly generated noise on top of the optimal 
sweep obtained from the power series, Eq. ©. The aver- 
age amplitude of the noise was varied between 0.01f^ a bi 
and lOOttjiabi- We then calculated the conversion effi- 
ciency of the power series detuning with the added noise. 
For noise amplitudes less than £ln a bi there was no notice- 
able change in the conversion efficiency or the dynamics 
of the atom-molecule conversion. While noise amplitudes 
larger the &,R a bi significantly effect the conversion dy- 
namics and efficiency. This indicates that as long as the 
size of the fluctuations are small enough to not shift the 
two states into or out of resonance with each other, there 
will be no observable effect. 




FIG. 1: (Color online) (a) Fraction of atoms converted into 
molecules as a function of log 10 r where r is the dimensionless 
magnetic field sweep time, r = T/TR a bi for fixed initial and 
final detunings A(-r/2) = -A(r/2) = -|A ] = ~50Q Ra bi- 
Here we chose Tn a bi = 1.6465 x 10 -5 s. The circles and squares 
indicate the conversion rate for the linear and optimal sweeps, 
respectively, (b) Landau-Zener factors, Slz, for linear sweep. 
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FIG. 2: (Color online) (a) The optimal sweep for r — 2. The 
dotted line is the optimal sweep obtained from Eq. © while 
the solid line is that obtained from Eq. (J7J. (b) Closeup of 
the optimal sweep obtained from Eq. JSJ for the power series. 



In order to understand why the optimal sweep is able 
to achieve such high conversion efficiencies even for very 
small t, we plot in Fig. 3 the instantaneous Landau- 
Zener factor S LZ {t) = fl 2 Rabi {t)/4\A(t)\ where n Rabl (t) = 
X|<^i(£)| and A(t) is the GA optimal sweep. At the begin- 
ning and end of the sweep 8hz{t) -C 1, which indicates 
that the sweep is non-adiabatic. However, at these times 
the detuning is so large that the atoms and molecules are 
decoupled. It is only for A(t) < ^R a bi(t) that the atom- 
molecule system is resonant and atoms can be converted 
into molecules. At these times 5jjz(t) > 1, which implies 
that the system is actually adiabatic . 

As one can see from Fig. 2, the optimal sweep rapidly 
goes from far off resonance to close to resonance and 
once near resonance the sweep rate slows down to main- 
tain adiabaticity. By using the narrow time interval at 
the beginning and end of the sweep where the detuning 
function changes the fastest, At, one can estimate the 
frequency bandwidth of the sweep, AE/fi rs 1/At, to be 
~ 4.5MHz. This is several orders of magnitude smaller 
than the typical energy spacing between molecular vi- 
brational levels in the closed channel. Consequently, the 
coupling between atoms and more than one molecular 
state can be ignored. 

At the same time the gradient of the detuning is large 
enough to prevent nonadiabatic oscillations, as one can 
see from Fig. 4(a). On the other hand the linear sweep 
results in nonadiabatic oscillations close to resonance 
that are the result of the coupling between the dressed 
states that is proportional to A(t)/Qn a bi- In the opti- 
mal sweep the atom-molecule conversion occurs gradu- 
ally over almost the entire sweep, while the conversion is 
limited to a very small window of approximately t rs to 
{« O.It for the linear sweep. The optimal sweep is sim- 
ilar to the switching scheme demonstrated in |2fj where 
the magnetic field was held at a fixed value close to the 
resonance for a certain time (t > 15ms). The switching 
scheme produced more molecules than the linear sweep 
because the atoms spent more time within the line width 
of the resonance where molecules can be created. How- 
ever, we emphasize that their hold times near resonance 
were several orders of magnitude longer than the nona- 
diabatic sweep times we consider here. 

To see whether the form of the optimal detuning was 
unique to Eqs. J3J and JSJ or were universal to two- 
state systems, we tested our GA on a two-level quantum 
mechanical system with the Hamiltonian, 



H = h 
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Fig. 5 shows the populations of the two levels \a) and \m) 
with all of the population initially in \a). The dynamics 
of the conversion of \a) into \m) is qualitatively the same 
as the atoms and molecules as is the shape of the optimal 
detuning function. This indicates that the general form 
of the optimal A(t) is generic to any two state system. 

The mean field energy shifts of the atoms (£%|0i| 2 , 
and Ui2\4>2\ 2 ) and molecules (C/22 102 1 2 and t r i 2 |</'i| 2 ) 
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FIG. 3: (Color online) Instantaneous Landau-Zener factor, 
<5iz(t) for the optimal sweep (solid line) and linear sweep 
(dot-dashed line). The sweep time is t = 2Tn a u 
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FIG. 4: (Color online) (a) Atomic, 



2 /no, (solid line) and 



molecular, |^| /no, (dashed line) fractions from the optimal 
sweep for r = 2. (b)Atomic and molecular fraction obtained 
from linear sweep for r = 1, 2, 6, 10. 



modifies the energies of the atom and molecules so that 
the difference in the chemical potentials between the 
atoms and molecules, A^i = ^molecule ~ Hatom, to zeroth 
order in \ is |3l| . 

A»=-hA(t) + U 22 \4> 2 \ 2 -U 11 \<f> 1 \ 2 + U 12 (\<p 1 \ 2 -\<j> 2 \ 2 ) 

(9) 
Consequently, even when A(t) = the energy difference 
between the atoms and molecules is not zero. 

In order to test the effect of the mean field shifts on 
the conversion process, we consider a range of possible 
values for the coupling constants [Zy to see how sensitive 
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FIG. 5: (Color online) (a) Populations of two level quantum 
system, a) (solid line) and \m) (dashed line) obtained from 
optimal sweep through resonance for r = 2. Also shown are 
atomic, |0i| 2 /no, (dotted) and molecular, 2\4>2\ 2 /no, (dashed 
dot) fractions for with Uu — U\2 = f/22 = and initial 
Rabi frequency chosen to match that of the two level system. 
(b)Optimal sweep, A(i), obtained from Eq. [S]for two level 
system (solid line). For comparison, the optimal detuning for 
the atom-molecule system is also shown (dashed line). 



the shape of the optimal sweep and the conversion effi- 
ciency is to variations of these parameters. We took the 
coupling constants to be either 0.1, 1, or 10 times bigger 
than the values given in Sec. II, which we denote here by 
\J\- . The maximum conversion efficiency for each of the 
27 different combination of values are shown in Fig. 6. 
The numbers in parentheses represent the relative mag- 
nitudes of the three coupling constants. For example, 
(1,10,0.1) means that U n = ll[°\ U 12 = 10U[°\ and 
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(0) 



22 



U 2 

In all cases, the optimal sweep shows very good con- 
version (> 98%) compared to less than 5% for the linear 
sweep even for very short sweep times. For instance, with 
Uu = lOt/f-L , the optimal sweep still shows better con- 
version than the linear sweep giving for example ~ 99% 
conversion compared to < 5% for r = 2. To understand 
this we note that from Fig. 4(a), most of the conversion 
takes place at the beginning of the sweep, t/r rs —1.0 to 
—0.5 where the atomic density is much larger than the 
molecular density. At these times the molecular energy 
lies above that of the atoms, —hA(i) > and the atomic 



mean field Uu 



shift reduces Afj,. 



As the Uij are varied, the qualitative form of the opti- 
mal detuning are the same as Fig. 2. There are, however, 
quantitative variations in the slope for different Uij as can 
be seen in Fig. 7. We tested the optimal detuning for 



FIG. 6: Molecular fraction as a function of non-linear inter- 
action terms for r = 2. The numbers in parentheses indicate 
relative magnitude of the mean field coupling constants Uu, 
U12, and U22, respectively. For instance, (10, 10, 0.1) means 
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and U22 = 0.1(7. 
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each of the the 27 combination of parameters denoted 



by A 



(J,- 



I (t) , on the other 26 combination of param- 



eters. In all cases the A 



(l,m,n) 



performed significantly 



better when applied to the other combination of param- 
eters than the linear sweep. These results imply that the 
optimal detuning calculated for any set of values of the 
two-body interaction can be used to convert atoms into 
molecules with greater efficiency than a linear sweep. 



V. DISCUSSION 

Very close to resonance, mean field theory breaks down 
due to the divergence of the scattering length and higher 
order quantum correlations can become important |28| . 
These higher order correlations have been shown to be 
important for time dependent experiments such as in 
Ref. j| where the atoms and molecules are placed in 
a coherent superposition state and allowed to evolve 
freely 32] . However, theoretical calculations that com- 
pare mean field theory and microscopic quantum calcula- 
tions that account for higher order correlations have been 
shown to give nearly identical results for the molecular 
conversion efficiency for downward magnetic field sweeps 
through a resonance |22,|23. Although these calculations 
were for linear sweeps, the agreement was very good for a 
range of A that included non-adiabatic sweeps for which 
Landau-Zener theory gives P mo iecules <C 1. These results 
support our use of mean field theory to model the atom- 
molecule conversion. It would, however, be of great in- 
terest to extend the work of Ref. J23 and compare mean 
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FIG. 7: (Color online) Optimal detuning function for several 
different values of Un, U12, and 1/22- 



field theory to quantum calculations for nonlinear mag- 
netic field sweeps. 

Throughout this paper we have ignored the effect of 
the decay of the molecules. If we assume a two-body in- 
elastic decay rate of 5 x 10 _11 cm 3 /s [Rj for molecule- 
molecule and atom-molecule collisions, we get a life- 
time of ti oss = 50/xs for the initial atomic density of 



4 x 10 14 e 



Assuming a simple exponential decay of 



the molecule number, we can estimate the fraction of 
molecules lost during the sweep by exp(— T/r; oss ). This 
estimate actually overestimates the molecule loss since 
it does not take into account the fact that t; oss is not 
constant but instead increases as the density of the gas 
decreases with time. exp(— T/ti oss ) should, neverthe- 
less, provide a lower bound on the fraction of molecules 
lost during a fast sweep. For T = Tr (,j = 16/is, the 
molecular conversion efficiency is reduced by a factor 
of exp(— Tjiabi/Tioss) — 0.73. Therefore, for a sweep 
time of TR a M the conversion efficiency using our opti- 
mal sweep and taking into account losses will be approx- 
imately eM-TR.au/Tioss) x 98.48% = 72%. 

In conclusion, we have applied a genetic algorithm to 
study the optimal time dependence of an applied mag- 
netic field, B(t), to convert atoms into molecules using 
a Feshbach resonance. Instead of the conventional linear 
sweep, we propose a nonlinear sweep through resonance 
that is characterized by a rapid approach of B(t) towards 
the resonance position followed by a gradual change in 
B(t) close to resonance and finally a rapid change in B(t) 
at the end of the sweep away from the resonance. This 
nonlinear sweep can result in almost 99% conversion of 
atoms into molecules even for sweep times as short as one 
Rabi period. We have shown that the qualitative form 
of this nonlinear sweep is independent of the specific na- 
ture of the two state system. We plan to extend this 
work by using a genetic algorithm to optimize molecule 
production via two-photon Raman photoassociation. In 
that case both the laser frequencies and intensities of the 
two lasers can be treated as time dependent quantities to 
be optimized. 
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